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SHAKE TEST OF RCTOR TEST APPARATUS IN THE 
40- BY 80-FT WIND TUNNEL 

Wayne Johnson* and James C. Biggers** 

Ames Research Center 
Moffett Field , California 


SUMMARY 

A shake test was conducted to determine the dynamic characteristics 
of a Rotor Test Apparatus on two strut systems in the Ames 40- by RO-ft 
wind tunnel. The rotor-off hub transfer function (acceleration per unit 
force as a function of frequency) was measured in the longitudinal and 
lateral directions, using a combination of broadband and discrete frequency 
excitation techniques. The dynamic data is summarized for the configurations 
tested, giving the following properties for each mode identified: the natural 

frequency, the hub response at resonance, and the fixed system damping. The 
complete transfer functions are presented, and the detailed test results 
are Included as an appendix. Finally, the report discusses the data 
analysis techniques developed to obtain on-line measurements of the system 
modal properties, including the damping coefficient and the damping ratio. 

INTRODUCTION 

A shake test was conducted to establish the dynamic characteristics 
of a Rotor Test Apparatus (RTA) in the Ames 40- by RO-ft wind tunnel 
(figure l). Of interest were potential resonances at the l/rev and 4/rev 
frequencies of rotors likely to be tested on the RTA, and potential ground 
resonance instabilities. 

The shake test was performed on the RTA module, without a rotor, on 
tw' strut systems in the wind tunnel, to determine the principle frequencies 
and damping of the structure. The rotor-off hub transfer function was 

*Research Scientist, U.S. Army Air Mobility Research and Development Laboratory 
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measured in the longitudinal and lateral directions : longitudinal, 

inplane acceleration of the hub due to longitudinal, inplane force; and 
lateral, inplane acceleration of the hub due to lateral, inplane force. 

With the hub transfer functions it ii possible to evaluate potential ground 
resonance and vibration problems of rotors to be tested on the RTA. 

The frequency ranges of interest are: 0-5 Hz for ground resonance, 3-7 «5 Hz 
for l/rev vibration, and 12-30 Hz for 4/rev vibration (based on a rotor 
speed range of 180-450 rpm) . The information required for each mode of the 
system is the natural frequency and the amplitude of the hub response, 
and for potential ground resonance modes we must know the fixed system 
damping as well. 


SYSTEM 

The system tested consisted of the RTA module, without a rotor, on 
the struts and balance frame in the 40- by R0-ft wind tunnel. The RTA 
module included the rotor hub, with the transmission locked, and two 
1500-HP electric motors installed (one of the motors was replaced by a 
dummy weight for this test). The total module weight was 30400 lb. 

Two strut/tip configurations were tested: a short strut system 

(8-ft struts with 5-ft tips) and a long strut system (15-ft struts with 
6-in tips). The short struts gave softer support of the module because 
of the flexibility of the tips. In the basic configuration the balance 
was free, with the scale system operating. The shake tests were also 
conducted with the balance locked, in order to obtain the cantilever strut 
modes. Finally, the system was tested with strut dampers installed, 
consisting of an extensible strut from the top of each main strut down to 
the rear of the balance frame, with a total of eight automotive shock absorbers 
as dampers. 


TEST APPARATUS 

A hydraulic shaker was attached to the blade grip of the rotor hub 
to excite the module by application of an inplane force, in the longitudinal 
or lateral direction. The other end of the shaker was attached to an 
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11600 lb reaction .ss suspended from a crane. Figure 2 shows the shake test 
configuration, for, lateral excitation. For longitudinal excitation, the 
shaker was attached to the forward blade grip, with the reaction mass 
over the module nose. The shaker servo control was operated in a stroke 
feedback mode. 

A load cell between the shaker and hub measured the applied force. 
Accelerometers on the hub measured the longitudinal and lateral response. 

Data were recorded for other accelerometers on the module and balance frame , 
but only the results for the hub response are presented in this report. 

The applied force and resulting hub acceleration data were analysed 
on-line to determine the dynamic characteristics of the system, using 
the Dynamic Analysis System (DAS, shown in figure 3)* The DAS is basically 
a time series analyzer and computer, utilizing Fast Fourier Transform 
techniques and associated software, and programs specific to this shake test. 

TEST PROCEDURE 

The frequency ranges investigated were 0-9 Hz for ground resonance 
and l/rev vibration modes, and 0-35 Hz for N/rev vibration modes. Broadband 
random input to the shaker was used, with a bandpass filter to shape the 
input spectrum. The low cutoff frequency was set at 0.5 Hz to avoid excitation 
of the reaction mass pendulum modes, and the high cutoff frequency was set at 
either 9 or 35 Hz to restrict the energy input to the frequency range of 
interest. 

The basic test plan, for each strut/module/ba lance configuration, 
excited in the longitudinal and lateral directions, was as follows. 

1) Random excitation, bandwidth .5-9 Hz; nominal force 
amplitude *200 and *400 lb (*400 lb point usually 
repeated ) . 

2) Random excitation, bandwidth .5-35 Hz; nominal force 
amplitude *100 and *200 lb. 
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3) Sinusoidal (discrete frequency) excitation at various 
force levels, at the low frequency resonances identified 
in #1» usually points were taken at frequencies near the 
resonance as well. 

There was some variation between runs of course. For future work, the use 
of narrow-band random excitation at each resonance would seen preferable 
to a sequence of (nominally) discrete frequency points as was the practice 
in this test. Narrow-band excitation offers the possibility of obtaining 
the data over the entire frequency range near the resonance in a single 
measurement. The reason for the narrow-band excitation is to concentrate 
the input energy into a particular mode, so for the highest force levels 
it may be necessary to narrow it down to essentially discrete excitation 
again. Still, the data from this test indicate that an accurate estimate 
of the system damping may be obtained from the single frequency point, even 
if it is not quite at the resonant peak (see the discussions below). 

The following six configurations were tested, with longitudinal 
and lateral excitation for each: 

1 ) Short s truts . 

2 ) Short struts, balance locked. 

3) Short struts, with strut dampers (B shocks). 

4) Long struts. 

5) long struts, balance locked. 

6) Long struts, with strut dampers (B shocks). 

ANALYSIS 

The data for the force applied to the hub and the resulting hub 
acceleration were analyzed on-line, utilyzing the DAS. The input signal 
f (force) and output signal a (acceleration) were sampled (digitized) 
at rate r, taking a total of N samples. The discrete Fourier transforms 
of f and a were calculated, and converted to engineering units using input 
conversion factors (lb/volt and g/volt). The products of the transforms 
gave the cross spectrum S^ o = F*A, and the input autospectrum S^ = 

S^ q and were averaged over K data records. Finally the transfer function 
of the hub response was calculated, from H = acceleration/force = 
averaged S^/averaged S^. 
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The computer searched the magnitude of the transfer function for 

resonant peaks. Then It calculated and printed for each peak the following 

quantities t the resonant frequency 0-3 ( Hz ) ; the magnitude of the force 

and acceleration at that frequency; the magnitude of the hub response |H| 

(g/lOOO lb and ln/lOOO lb); the phase of the response, /»H (deg); the 

fixed system damping coefficient C (lb/fp6), calculated from H at the 

8 

resonant frequency; the damping coefficient, modal mass, and damping 
ratio (C g , M, and ^ ), calculated from integrals of H through the resonant 
peak; and the damping ratio £ , calculated by a least-squared-error 
parameter identification technique from the data for H near the peak. 

In addition, ground resonance parameters (critical rotor speed and required 
lag damping) were calculated, for a particular rotor. 

The magnitude of the transfer function, I Hi vs • V) , was displayed 
on a CRT. A picture was taken as a record of the complete transfer function. 

The discrete frequency excitation points were analyzed in the sane 
manner. However, the response was only evaluated at the single line corresponding 
to the input frequency. 

Further details of the analysis techniques are given in appendices t 
a discussion of the discrete Fourier transform (Appendix A); the local 
maximum discriminator (Appendix B); calculation of the fixed system damping 
from the transfer function (Appendix C); LSE parameter Identification of 
the dam ling ratio (Appendix &) ; and calculation of the damping from integrals 
of the transfer function (Appendix E). 

The following parameters were used for the analysis of the data 
in this test i 

l) Ground resonance and l/rev dynamics (0-9 Hz)« sample rate 
r * 20.48/sec, number of samples N ** 512, number of records 
K ■ 10; total sample time T * 250 sec, spectrum frequency 
increment ■ .04 Hz. 
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2) N/rev dynamics (0-35 Hz); sample rate r - 81.92/sec, 
number of samples N * 256 , number of records K * 20; 
total sample time T • 62.5 sec, and spectrum frequency 
increment Asa * .32 Hz. 

For future work, it would probably be better to take more records for the 
35 Hz bandwidth excitation, to further reduce the noise in the data; 

K * 40 (hence T • 125 sec) should be about right. The use of Hanning to 
smooth the data is usually recommended (see the references of Appendix A), 
but it was only occasionally used in this test. 

RESULTS 

The results of this test are the dynamic characteristics of the six 
configurations investigated, specifically , the frequencies and response 
amplitudes of the principal modes identifiable in the hub transfer functions. 

Figure 4 demonstrates the repeatibility of the transfer function 
measurements. It shows three separate measurements of the longitudinal 
response on the short struts. There is excellent correlation between 
the three points. Figures 5 through 10 present the transfer functions 
for the six configurations tested. The lateral and longitudinal hub 
responses are shown, in the 9 and 35 Hz excitation ranges for each. 

The abscissas in the figures are frequency, from 0 to 10 or 50 Hz, and 
the ordinates are the magnitude of the transfer function in g/1000 lb. 

Tables 1 and 2 summarize the dynamic characteristics of the six 
configurations tested. The tables give the following quantities for each 
of the longitudinal and lateral modes identified* the resonant frequency 
W (Hz); the magnitude of the hub response H (g/1000 lb and in/1000 lb ) 5 
and, for the potential ground resonance modes, the fixed system damping 
coefficient C g (lb/fps). The hub response and damping coefficient data for 
the long struts, lateral shake, balance free and balance locked (runs 1 ? and 
18) are somewhat uncertain because of a problem with the accelerometer 
calibration. However, the conversion factor (g/volt) used for these two 
runs was certainly within 23& of the correct factor. The frequency and 
damping ratio data are nof affected by this problem. 
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TABLE 1. Summary of Dynamic Characteristics t Short Struts 
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In the lateral response on the short struts (balance free, no 
dampers) we observe two close modes at the lower resonance, around 2 Hz. 

These are the balance side and balance yaw modes, involving considerable 
module yaw and side motion as well for this case. Figure 11 shows the details 
of the two modes, expanding the magnitude and phase of the transfer function 
in the range 1-3 Hz. The detailed test data are given in Appendix G, 

Table Gl. The dynamic situation is as follows. For the short strut 
configuration, the uncoupled balance lateral modes and cantilever strut 
side mode have about the same frequency, around 2.4 Hz (see Tables 1 and 2). 
Thus there is considerable coupling of the balance and module motion for 
the complete system, with the typical behavior that the frequencies of 
the coupled modes are driven apart. The balance mode frequencies are 
decreased, the strut mode frequency is increased, and the damping for the 
balance modes is reduced. For the long strut configuration, the lateral 
cantilever strut mode frequency is around 3*5 Hz, well above the balance 
mode frequencies (see Table 2). Therefore the balance mode frequencies 
are not lowered significantly for this configuration (note that the balance 
side mode frequency is expected to be about times the balance 

longitudinal mode frequency, since there are two side force scales and one 
drag force scale in the balance system). 

The test data show a nonlinear behavior for the damping of the balance 
modes. The damping for high excitation level and high response amplitude 
consistently was significantly lower than the damping measured at low levels 
(the data in Tables 1 and 2 are the values for low excitation level). 

Figure 12 shows the general trend, fo~ all the configurations tested. The 
ratio of the balance mode damping to its value at low excitation levels 
correlates well with the rms value of the exciting force. The linear 
range extends up to 30 or 40 lb (rms), and at high excitation the damping 
levels off at about 40 to 50# of the low excitation value. The detailed 
test data are given in Appendix G. A correlation between the excitation 
level (broadband and discrete) and the balance drag scale motion for the 
longitudinal balance mode is also given in Appendix G (run 11, Table Gl). 
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From the frequencies of the modes we may identify l/rev and N/rev 
resonances for the operating range of a particular rotor, and with the data 
on the magnitude of the hub response assess the vibration potential of 
these modes. From the frequency and fixed system damping we may assess 
tho ground resonance stability of articulated and soft-inplane hingeless 
rotors on this Rotor Test Apparatus. A simple ground resonance stability 
criterion, giving the critical rpm ranges and the lag damping required 
for stability, is discussed in Appendix F. More detailed calculations 
of the dynamic stability are recommended however. 

The tables of Appendix G present in detail the shake test dc ta 
for the six configurations investigated. 

DISCUSSION OF ANAL f SIS TECHNIQUES 

Several methods were used to calculate the modal parameters from 

the measured transfer function. The quantities required are i t.ie 

natural frequency , damping coefficient C , damping ratio t , and 

s ^ 

modal mass M (note that these parameters are related by C g - 2 J*V <)• 

The natural frequency was estimated using three points around the 
experimental peak (Appendix D) . This technique gave satisfactory results. 

The damping coefficient was calculated from the transfer function 
at a single frequency point (Appendix C), and from integrals of the transfer 
function through the peak (Appendix Ej this method was used only for runs 
17 and 10, Tables G4 and G5). Both methods worked well, and the two teckmiques 
gave comparable estimates. The experimental data (Appendix G) for the single- 
point calculation of C^ during discrete frequency sweeps near a resonance 
demonstrate that this method gives an estimate of the damping which is 
indeed relatively insensitive to frequency, l.e. roughly constant in the 
vicinity of each peak (see Appendix C). The integral method of calculating 

C is less sensitive to noise in the transfer function data, but for very 

s 

close modes one must watch that the limits of integration cover only one 
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resonant peak. With discrete or very narrow-band excitation, only the single- 
point estimate of C is applicable of course. The use of both methods is 
recommended to obtained the best extimate of the damping coefficient. 

To calculate the damping ratio and modal mass (£ and M), the LSE 

parameter identification techniques described in Appendix D were used , 

with four iterations after the initial estimate of the parameters, and either 

5 or 10 points for the curve fit '.round the resonant peak. For an ideal 

transfer function (no noise in .e data) these techniques worked well, 

especially the two-parameter algorithms. For real data however, i.e. an 

experimental transfer function measurement including noise, the methods of 

Appendix D were not satisfactory The one-parameter algorithm ( £ from 

\*\ ) did no better than the initial estimate of £ from three points. 

The two-parameter algorithms ( and M from either )h| or H) either 

gave little improvement over the initial estimate, or simply did net 

converge. The difficulty is probably due to the fact -that the derivatives 

in the iteration formulas are singular at £ m 0. Theie is the possibility 

of better success using an algorithm to identify C and M from the transfer 

s 

function. 

The damping ratio and modal mass and M) were also calculated 
(for runs 1? and 18, Table G4 and G5) from integrals of the transfer 
function, as described in Appendix E. This technique worked well, and 
its continued use is recommended. 


- 11 - 



APPENDIX A 

The Discrete Fourier Transform 


1 . References 

Bendatf Julius., and Piersol, Allan G., Measurement and Analysis 
of Random Data . John Wiley 4 Sons, Inc., New York, 1966 

Jenkins, Gwilym M., and Watts, Donald G., Spectral Analysis and its 
Applications . Holden-Day, San Francisco, 19w 


2. Definition and Application 

The input signal (force, f) and output signal (acceleration, a) 
are sampled (digitized) at rate r, until the total number of samples N 
is collected. The result is a discrete time series of data, at t = n At, 
n = 0...N-1 ( A>t ■ l/r, with sampling period T = N/r). The discrete 
Fourier transforms of the input and output are calculated, using Fast 
Fourier Transform (FFT) techniques, according to the expressions 


X(0 — £L * e, 

^ n* o 




The result is a discrete spectrum, at the N/2 frequencies u = k Atd , 
k = 0...(N/2-l) ( AtO = r/N = l/T Hz, with a maximum frequency — spectrum 
bandwidth — oi ca = r/2 Hz). 

2HSLX 


The input and output transforms are multiplied then to obtain the 
cross-spectrum S iQ = F * A and the input autospectrum = F * F. The 
spectra and are averaged over a total of K records of data. Then 
the system transfer function is calculated ass 

H = a/f = average S ^/average 


3. Relation to Continuous Fourier Transform 

The Fourier transform of a continuous, non-periodic time function 
x(t) is defined as 


i 

S A I 


— t u*tr 



The discrete Fourier transform makes two approximations! first, a finite 
length record of data is transformed; and secondly, a finite number of 
samples are taken during the record. 

With a finite length record, it is assumed that the data is periodic 
outside the record; hence we calculate the Fourier transform of a periodic 
functions -j- • » 

S 0 * t0< ^'“ a*- 

at the discrete harmonics 2tCk/T. This may be considered the Fourier 

transform of x(t) times the window w(t) which is open only for t = 0 to T, so 

X,( w > — 2 ^- S *!*!'*• to *T" ' a* 

- bo 

(using the convolution theorem). The time window is 

2 Vt < *r 

o otherwise 

so the frequency window (the Fourier transform of w) is: 

V , osr/z tf»V2 

W = -ZSyT * 

which has amplitude 1 and bandwidth = 2*nr/T. Thus 

c-r* 

This Aft is the same as the frequency increment in the discrete spectrum. 

So each line in the discrete transform may be viewed as the integral of 
the continuous transform over the interval o — Aft/2 to <ft+Av»/2. . 

With only a finite number of samples in the record., we calculate 
as an approximation to the discrete transform! 
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£ *(** c 


i 2ir kA AX/r" 



At 


(the summation being the discrete approximation of the integral). Since 
At/T » 1/N, this is identical to X(k) defined in section 2 above. 

The finite length record means that only the discrete harmonics « k A** 
of the transform are calculated. The finite number of samples means that 
the maximum frequency of the spectrum is *" i 1Iiax “ N/2 * “ r/2 (the Nyquist 
frequency). It is necessary to filter the analog input and output signals 
with a low-pass cutoff frequency at or below i \ ax * in order to avoid aliasing 
of the discrete spectrum by harmonics above the Nyquist frequency, which 
can not be discerned by sampling at the discrete rate r. 


4. Noise 

Because of process and measurement noise, we do not calculate a 
deterministic spectrum, but rather a statistical estimator of the spectrum. 

In order to reduce the noise in the estimate of the spectrum, it is necessary 
to average the data. Thus we take K records, and calculate the average 
spectrum . 

s « t 


This sample spectrum has an unbiased mean, and a variance of 

Nour S. _ ~ 

* “ k 

The standard deviation is thus inversely proportional to K 2 (compare with the 
similar result for the standard deviation of a sample mean). The total 
sample time is KT = K/(AtO , so for a given time it is necessary to compromise 
between the accuracy of the data and the frequency increment in the spectrum. 
Bendat and Piersol suggest using a minimum of K = 10 records. The statistics 
of the transfer function H (the ratio of the average cross spectrum to the 
average input autospectrum) are more complex (the reader is directed to the 
references given above), but the K~ 2 behavior of the spectra is sufficient for 
the present purposes. 
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5. Choice of Parameters 

The parameters r, N, and K are required to define the sampling and 
averaging process in the analysis of the data. For a given bandwidth of 
the data, the sample rate r suggested is 

r = 2.5 * bandwidth data 

(*■> ■ r/2 ■ 1.25 * bandwidth). A low pass filter on the signal is also 

required, to avoid aliasing in the discrete transform. The number of 
samples N is then chosen from r and the required frequency increment in 
the spectrum AO , as N = r/A*-» (FFT routines used require also that H 
be a power of 2). We choose AtJ to define the resonant peaks sufficiently, 
from 5 ua '*V2. (which gives about 5 points covering the £ power 

bandwidth of the peak; is the natural frequency and £ the damping ratio 
of the mode). Finally, the number of records K is chosen for the desired 
accuracy (noise level) of the spectrum. At least 6 to 10 records are 
desired; the principle restriction of the number of records is the total 
sample time KN/r. 
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APPENDIX B 

Local Maximum Discriminator 


1 . Problem 

It is necessary to identify the resonant peaks (i.e. the natural 
frequencies) of the experimental transfer function. The experimental 
transfer function has measurement and process noise however , so it is not 
possible to identify the peaks by simply searching for all the local maxima 
of the data. An algorithm must be developed which will discriminate the 
true peaks from the spurious local maxima due to noise. 


2. The Algorithm 

We have the data for the magnitude of the transfer function, which 
may be written H g = H + h, where H is the true value and h is random noise 
in the measurement. Assume h has a normal distribution with zero mean and 
standard deviation v* , hence probability distributions 




■ns? 


1 

Assume = H/K 2 , where K is the number of records of data in the average 
of the cross spectrum and input autospectrum calculated to find H (see 
Appendix A). 


Consider then the probability of a peak at a certain frequency 

i.e. the probability that Hjj-H^O for all nearby frequencies. This is 

the probability that h h„ - AH , where AH = H - H ; which is : 

N e e ew e 






s 


ss 




<&JU 


bo 


J 
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This Integral nay be expressed in terns of the error function. 

The product of the probability Pr evaluated at several points around 
iOjj is the probability that all local values of H are less than the H at *v 
Therefore we take as a discriminator of the local naxima the parameter 

c - c •fr O* 

C is evaluated at all frequencies of the transfer function. If C is above 
a certain confidence level for any frequency, we consider that point a 
resonant peak of the transfer function. 

The parameter C has the following properties. For a local maximum , 

C 1 , while C is near 0 for a local minimum. If AH = 0 for all points (i.e. 
the experimental data constant), then C * Finally, with &\J*~ = 1, 2, or 3 
we obtain C = .76, .92, and .98. 

3. Application 

For on-line evaluation of the data (locating the resonant peaks of 
the transfer function and calculating the system properties there), it 
is better to use a rather low confidence level on the discriminator (so a few 
false peaks are located, which are easily discarded by the engineer), rather 
than to use a high confidence level which will occasionally miss a true 
peak because of excessive noise. It is also found that the parameter C 
is a more sensitive discriminator of the peaks if many points are used 
to evaluate C for each frequency. 

For the present test, a confidence level of 65 to 70 (C above the 
confidence level considered an indication of a resonant peak) was satisfactory. 
The parameter C was calculated using 12 points (**m" in the definition of C 
above) around each frequency. 
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APPENDIX C 

Fixed System Damping from Transfer Function 


To evaluate the ground resonance stability of a rotor on 
a flexible support, it is necessary to know the damping coefficient 
of the modes. This may be obtained from the hub impedance by the following 
method. Consider the mass/spri ng/damper systems Mx + C g x + Mc^x- f. 

The response of the hub acceleration to the applied force is the transfer 
function 


Ll. i - 


— 


(**« — w 1 ) -t- 


where is the natural frequency, M the generalized mass of the mode, 
and C g the damping coefficient. It follows that 


or 


5 ~ \Y\\ X 


c.s = ns -A 


p v\ 
\vu^ 


with the dimensions [y*J- Hz, [] H J - g/1000 lb, [c "3 = lb/fps. This 
is the expression used to calculate the damping of the rotor support, 
from the experimental measurement of the hub response. 


At the resonant frequency ( l» * ) this result becomes C ” *°/ \ h\ . 

s 

In general the previous form is preferable however, since it holds for all u , 

not just at the peak. Thus it is possible to evaluate C even though the 

s 

calculation is not performed exactly at the peak (for multimode systems it is 
necessary to be at least close to the peak of course). The experimental date*. 
(Appendix G) shows that the damping calculated by this expression is quite 
consistent in the vicinity of the resonance of each mode. 
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APPENDIX D 

Least Squared Error (LSE) Parameter Identification 
of Damping Ratio from Transfer Function 

1. LSE Parameter Identification 

It is desired to fit an analytic function H(W ,u^) — where t*> is 
the frequency, u^ are free parameters (e.g. the damping ratio ^ ), and 
H may be either a complex transfer function or the magnitude — to experimental 
data H .(<v at the discrete frequency points . We shall find the 
parameters u^ to minimize the squared error 

For complex H, this error is the sum of the distances between H and H on 

© 

the complex plane (RmHvs. Im H , i.e. the phase plane). The minimum € 
is given by the solution ofi 

= ± \ ® o 

k 


If H is linear in the parameters, the above is a set of linear algebraic 
equations which may be solved directly for the parameters u^. In the present 
case however, H is not a linear function of u^, so a solution try numerical 
methods is necessary; we shall use Newton's method. From 


£ 

dui 


^5 t 

a»Al \ ^ 




z 




'r 




->c*> 


it follows that the Iterative solution of 


^ — w 


t 


bf/& u t = 0 is 



where u is a vector of the parameters u^ (nth iteration), and the derivatives 
of f are evaluated using 
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, hence the solution of the parameter identification 


Here f - Z I H - H I 2 
y. 1 e ■ 

problem ist 


i* ar V* 





k 


2. Transfer function 

We shall fit the measured transfer function in the neighborhood 
of a resonant peak to the theoretical transfer function of a mass/spring/ 
damper system. Considering the acceleration response to an applied force, 

the transfer function is 

^ — u> x /wt 

H ^ ~ + i» 2^^* 

Note that in general the parameter m (mass) is a complex number, because 
it accounts for the influence of other modes of the system in the vicinity 
of any particular resonance. The magnitude of H is: 



Fitting H to the experimental data around a peak requires the identification 
of four parameters then: the damping ratio £ , the natural frequency Li*, , 

the mass |vt>\ , and the phase angle /.M (only the first three are involved 
in fitting the magnitude of H to the experimental data). 

Because of limitation of computer core and language, we consider only 
the Identification of one or two parameters. The following cases will be 
considered in detail: fitting |h^ to |H e | by identifying £ 5 fitting 

) H \ to | H | by identifying T and \u\ ; and fitting H to H by identifying 

£ and . An initial estimate of the parameters is required to start 

the iterative LSE solution. It is assumed that the initial estimate of the 
parameters not corrected by the LSE solution (in particular the natural 
frequency u^) is satisfactorily accurate. 
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where 

£w ss I Hc u /Htp | 

* I I 
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4 , Damping ratio from 1 H I 

The LSE iterative solution is 

* C.u + 

or 


S' ~ £U-UeV- 

J-fy c*-^y 




where 


r •= 


C jr^ CU-He^ 


£ 

K 


3 S 


5 




r W* b" Va ("We-K - ) 

__ 14 — 

% J“ Ci2c\Ot>' s/2 CU t - ±W) 

k J J 

b ~ ) 2 

, X 

Y\ = — b~ ** 

yv\ 

5. Damping ratio and Maas from f H 1 

With i/|ml and D = ( M l -uai ) 2 + (2£tJwu) 2 , we have 

H *=y> tO*D 2 . The derivatives required are « 

A = l IjCH-iO* = f CWe-»V wi b" 3/i f| 

B= £ fHe-H^ = C(U t -H > ) (~Zu l 
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6. Damping ratio and Mass from H 

Using the initial estimate of and 6S* , we shall match the 

experimental data H e e*^ m to the complex transfer function 

H - 7«w* b"' 1 C — 1 2 *) 

where D = (fc» — !•»* ) + (2£tf»d0 and y* * l/)m\ . Then the squared error is 

C s t ^ ( o< -c^u.'' ^ ^ 5 > tr' + 1*0*] 

where 

cm * 2 fccHt, w l 

2 w l (-2uO 
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The derivatives required ares 

» - £ AtlM-Hel* 

- ?/* a 6 *5 -j 

fc* e. s S 2w»<r' 

K. VC 

and the LSE iterative solution isi 

( £ \ _ / £ \ _ ! /AB-S»\ 

V /* "" V/e ' oJ3> — C Z V fco — Ac / 
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APPENDIX E 

Damping Ratio from Integral of Transfer Function 


The damping ratio, mass, and damping coefficient may be calculated 
from integrals of the system transfer function. This method is an alternative 
to the single point or curve fit techniques described above (Appendices 
C and D). Assuming a single mode transfer functions 

_ 2 . _ — 

~ 5 pl-wj*. + i Zjjuvo*. 


it may be shown that the damping coefficient and mass are given byi 


s: 


nr/z 

)o w 

and then the damping ratio is C - C /2mu) . 

J 8 R 

To apply this result to experimental data, the transfer function is 
integrated through each mode from 0.8u) n to 1.2u^« Correcting for the 
finite limits, we obtains 

1510 


5 


.33 





u> 
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e s = W5-'i 


with dimensions M - g/1000 lb, [«3 * Ha, [■] » lb, and [ q ^\ = lb/fps. 
Note that the result for C is independent of the limits of integration} 
compare with the expression in Appendix C. For extremely close modes it 
may be necessary to integrate over a smaller range around U n ; the above 
limits were satisfactory for the present test however. The natural frequency 
U> n may be obtained from the three-point curve fit around a local maximum, 
as described in Appendix D, part 3 » 

By calculating the system parameters from integrals of the transfer 
function, the effect of noise in the experimental data is reduced. 

However, the above expressions are not unbiased estimators of ^ and C g . 

With the factor |h|^ in the denominator, the calculation of £ and C g in 
the presence of noise will underestimate the true values. The error 
in the estimate will be of the order K , where K is the number of data 
records over which the spectra are averaged (see Appendix A). The estimate 
is conservative at least, and for the present cases the error is only 5 to 
10$. if desired, the calculations of £ and C g may be multiplied by 
(K+l)/K as an approximate correction for the bias error. 
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APPENDIX F 

Ground Resonance Stability Criterion 


1 . References 

Coleman, Robert P., and Feingold, Arnold M., "Theory of Self -Excited 
Mechanical Oscillations of Helicopter Rotors with Hinged Blades," 

NACA Rept. 1351, 1958 

Deutsch, M.L., "Ground Vibrations of Helicopters," Journal of the 
Aeronautical Sciences , vol, 13. no, 5, May 19^6 

2. Ground resonance 

Ground resonance is a mechanical instability involving the coupled 
dynamics of the rotor lag and hub inplane motion. An instability is possible 
at the resonance of the low frequency lag mode (frequency SZ ) and 

a fixed system mode (frequency lO* ) , if the product of the fixed system 
damping and the rotor blade lag damping is below a critical level dependent 
on the blade inertia and leg frequency. 


proximate stability criterion 
The critical rotor speed for resonance is 




and the system damping required at resonance is 



where to* = support natural frequency 

= support damping 
- blade lag frequency (rotating) 

= lag damping 
N = number of blades 

S- = first moment of blade mass about lag hinge 
5 (i.e. mass * radial C.G. location) 
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With dimensions £si3 = rpn, = Hz, “ P® r rev * C0 X "1 * lb/fps, 

C Cg*} = ft- lb/rad/sec, and [s^*3 * slug-ft, the stability criterion is» 








Usually the required lag damping for stability (C is increased by a 
margin of 30 to 50% to obtain an engineering estimate of the stability 
boundary. 


This criterion is based on the assumption of a small ratio of blade 
mass to the rotor support mass, which is visually quite true. For extremely 
small fixed system damping this approximate criterion may not be conservative 
however. In general a detailed analysis of ground resonance stability is 
recommended. 
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APPENDIX G 

Rotor Test Apparatus Shake Test Data 

The tables of this appendix present the data for the resonant 
frequencies of the hub transfer functions (lateral acceleration due to 
lateral force, and longitudinal acceleration due to longitudinal force). 

The following configurations were tested: 

Table Gl. Short struts 
Table G2. Short struts, balance locked 
Table G3* Short struts, with strut dampers (8 shocks) 
Table G4. Long struts 
Table G5- Long struts, balance locked 
Table G6. Long struts, with strut dampers (8 shocks) 
The following quantities are given in the tables: the resonant frequency 

40 (Hz); the amplitude of the hub response H (g/1000 lb and in/1000 lb); 
the phase of the response ZH (degrees) ; the fixed system damping of the 
mode C (lb/fps); the damping ratio t (per-cent of critical damping); and 

5 

the amplitude of the exciting force F (rms lb, with "D w indicating 
discrete frequency excitation)* Several sweeps of discrete frequency 
excitation in the vicinity of resonances were made, and the data are given 
for the entire sweep as well as for the peak. 

The hub response and damping coefficient data for the long struts, 
lateral shake, balance free and locked (runs 17 and 18, Tables G4 and G5) 
are somewhat uncertain because of a problem with the accelerometer calibration. 
However, the conversion factor (g/volt) used for these two runs was certainly 
within 2% of the correct factor. The frequency and damping ratio data 
are not affected by this problem. 
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Table G1 . (continued) 
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Table Gl* (continued) 
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Table G2, Short struts, balance locked . 
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Table G2. (concluded) 
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Table G3. (continued) 
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Table G3 * ( continued ) 
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Pable G3 . ( concluded ) 
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Table G4. (continued) 
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Table G^. (continued) 
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Table G4. (continued) 
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Table G?. Long struts, balance locked 
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Table G5. (concluded) 
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Table G6. (continued) 
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Table G6. (continued) 
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Table G6. (continued) 
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Table G6. (concluded) 
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Figure 4. Fepeatibillty of transfer function measurement. 
Short struts, .5-9 Hz broadband excitation, 
longitudinal hub response. 
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»5-9 Hz broadband 
excitation, longitudinal 
hub response 

(run 10, point 13) 
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hub response 

(run 10, point ?) 



Figure 6. (concluded). 
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hub response 
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.5-35 Hz broadband 
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Figure 10, Long struts, with strut dampers (H shocks). 
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Figure 12. Nonlinear damping characteristics? reduction 
of fixed system damping (C ) with excitation 
force level (F), for modes involving -oredo.ninantly 
balance motion. 
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